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Abstract 

We study the different nucleus-nucleus collisions, O+Au, S+S, S+Ag, S+Au 
and Pb+Pb, at the CERN-SPS energy in a one-fluid hydrodynamical ap- 
proach using a parametrization based on baryon stopping in terms of the 
thickness of colliding nuclei. Good agreement with measured particle spectra 
is achieved. We deduce the mass number scaling behaviour of the initial en- 
ergy density. We find that the equilibration time is nearly independent of the 
size of the colliding nuclei. 
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I. INTRODUCTION 



From heavy ion collisions at ultra-relativistic energies one would like to learn, e.g. how 
high energy and baryon number densities can be reached in these collisions and to what 
degree the produced matter is equilibrated. Answers to these questions would tell us which 
regions in the nuclear phase diagram are accessible in nuclear collisions. It would also show if, 
and to what extent the quark-gluon phase has been reached in present SPS experiments [0. 

There are no direct measurements of the achieved energy and baryon number densities 
and therefore one uses models to extract the information from the experimental data. One 
commonly used approach is hydrodynamics || which correlates the initial densities with 
measured spectra of final particles. By definition, hydrodynamics assumes local thermal 
equilibrium in the expanding matter and is therefore unable to answer the question of 
thermal equilibration. 

Our main goal in this work is to investigate the scaling of initial energy density in heavy 
ion collisions as a function of the size of the colliding nuclei at (nearly) fixed energy. Therefore 
we study and compare five different collisions at CERN-SPS energies, O+Au, S+S, S+Ag, 
S+Au and Pb+Pb, using the same parametrization for the initial state. In Section || we 
introduce the parametrization of the initial conditions. Section |T| discusses the results and 
in Section |TV| conclusions are drawn. 



II. INITIAL CONDITIONS 

At ultra-relativistic energies nucleus-nucleus collisions are expected to become transpar- 
ent, meaning that the matter of final particles is produced in a state of initial collective 
motion. We have not succeeded to reproduce, even qualitatively, the baryon rapidity distri- 
butions if we assume that the produced matter is initially completely stopped, commonly 
referred to as the Landau initial conditions. As an example we show in Fig. p] the net-proton 
distribution in a lead-on-lead collision as obtained by assuming that the matter of the initial 
fireball is at rest. The calculation is compared with the measurement of NA49 ||. The 
disagreement between the data and the calculation is, as expected, larger for lighter nuclei. 
This indicates that the primary collision stage cannot be treated within one fluid hydrody- 
namics since this would lead to compressed matter at rest. Instead, we start the simulation 
after the system has equilibrated [fU|. To use this procedure we need a parametrization of 
the initial conditions. Little is known about the initial state and a chosen parametrization 
can be justified only by comparing with the data. The parameter space of the initial condi- 
tions, which give reasonable agreement with data turns out to be large. Partly this results 
from the freedom in choosing the equation of state. Therefore it is difficult to pin down, 
e.g. the initial temperature achieved in a heavy ion collision. 

Here we exploit the idea to fix the initial conditions using the local thickness of colliding 
nuclei. We parametrize the nuclear stopping locally in the transverse direction in terms of 
the nuclear thickness function and then constrain the energy density from conservation of 
energy. The advantage of this approach, when combined with the geometry of the colliding 
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nuclei, is that fixing of the parameters fixes the initial conditions for all different collisions 
at the given energy. It also gives a better basis to compare different collisions and to draw 
conclusions on the dependence of the initial state on the nuclear size. 

The main features of our parametrization of the initial conditions are described in ||. 
Here we summarize the most important ingredients. Note that we consider only impact 
parameter zero collisions i.e. azimuthal symmetry is assumed. 

One of the difficulties in hydrodynamic simulations is the treatment of the surfaces. 
E.g. in the transverse direction the outer edges never behave hydrodynamically and even in 
head-on collisions a certain fraction of nucleons never suffer a single collision. These nucleons 
are not detected in the experiment and their energy does not contribute to the production 
of final matter. To account for this we use effective nuclear sizes, i.e. we replace the mass 
numbers of incoming nuclei, A and B by 

A eS = £A , B cS = £B (1) 

and fix the geometry in terms of these effective mass numbers A eS and B cS . In fixing the 
value of £, always close to one, we also compensate for the few percent losses of baryon 
number and energy in our numerical code. 

We start out by parametrizing the rapidity distribution of baryons dN^/dy as a sum of 
projectile and target contribution for symmetric collisions ||, 



dy dy dy 

C p exp(ax y + bx y + cx y ) + C T exp(—ax y + bx y — cx 



p 



using a variable x y = y/y ma , x , the rapidity scaled with y max = y^ m , the projectile rapidity in 
the cm frame. The factor (1— x y ) ensures that the distribution goes to zero at the boundary 
of the phase space when effects like Fermi motion are neglected. The functional form is 
motivated by the experimental proton rapidity distributions in p + p collisions || and the 
possibility to control the amount of stopping. C p and C T are normalization constants. The 
main idea of our parametrization is that the above initial distribution is determined locally 
in transverse plane. This is implemented by defining the parameters cl(Ta), b(T^) and c(Ta) 
as functions of the local nuclear thickness 

t a(p) = J dz n B (z, p) , (3) 

where n&(r) is the nuclear density for a nucleus of mass number A, z is the longitudinal 
and p the transverse variable: f — (z,p). We use the Woods-Saxon parametrization for the 
nuclear density 

nB (f) = exp[(|r1-iL)/a R ] + l ' (4) 

with 

R A = 1.12 fm x A 1/3 -0.86 fm x A" 1/3 , (5) 
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clr = 0.54 fm and no = 0.17 fm 3 0. The functional dependence of a, b and c on T A is 
chosen to be || 

a{T A ) = 1.5 {(JppTA)- 1 
b(T A ) = f3 a (1 - a P pT A ) 

c{T A ) = 3.0 , (6) 

where a pp = 32 mb is the total inelastic cross section for p + p collisions at SPS energy. The 
first motivation for this parametrization is its simplicity. Second, in the case of one collision, 
i.e. (J pp T A = 1, the experimental proton rapidity distribution for p + p collisions is recovered 
|J. For a pp T A > 1, the b coefficient is negative indicating increasing stopping with growing 
nuclear thickness, the strength being controlled by f3 s . It is fixed from fits to baryon rapidity 
spectra in heavy ion collisions and therefore its determination includes the hydrodynamical 
evolution. Changing the hydrodynamical evolution by choosing different initial energy and 
(longitudinal) velocity distributions or equation of state results in different optimal values 
for j3 s . In the analysis presented here the value (3 S = 2.25 adequately reproduces all results. 

We would like to emphasize that T A (p) depends on the transverse coordinate p and 
therefore stopping at the center is different from stopping at the edge. As an example we 
plot in Fig. H] the baryon density 777,3 in the cm frame for S+S collision. One sees clearly the 
increase in transparency with radius. This makes the one-fluid hydrodynamical treatment 
of nucleus-nucleus collisions more realistic. 

For the initial local energy distribution we choose a Gaussian form 



■{y-yof 

2al 



(y/Vmax) , (7) 



dy 

where the width a e and the normalization C £ depend on the transverse coordinate p. The 
normalization C £ is given by energy conservation and the value yo is identified with the 
center-of-mass rapidity of the collision locally in the transverse plane. Since our hydrody- 
namic calculation is performed in the overall center-of-mass frame of participating nucleons, 
yo(p) = for zero impact parameter collisions of equal nuclei. For asymmetric collisions the 
thickness of target and projectile varies differently with transverse radius p and therefore 
the rapidity of the local cm frame yo(p) changes in the transverse plane. 

In Eq. (0) the only free parameter is the width <y e which is taken to depend on the nuclear 
thicknesses as [|J 

° AP) = [a pp T A (p)al p T B (p)} * ' (8) 

where a pp is included to make the denominator dimensionless. The constants c e and a e are 
determined from an overall fit to the investigated collision systems with the result c £ = 0.8, 
a £ = 0.125. Q 



1 These values are slightly different from those used before in Q due to the extension to more 
collision systems. 
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So far baryon density and energy density distributions are specified only as functions of 
velocities. To convert them into spatial distributions we have to specify the velocity field. In 
the case of the Bjorken model 0, the scaling ansatz for the longitudinal velocity is v z = z/t 
and initial conditions are usually defined at fixed proper time To- Since scaling cannot hold in 
a finite system, we have, for numerical convenience, chosen the rapidity y instead of velocity 
v z to have the linear z-dependence: 

y(p, z) = k{p)z , v z (p, z) = tanh[/«(p)z] , (9) 

where the proportionality constant k depends on the transverse radius p ||. For small z 
this ansatz can approach the velocity profile of the scaling solution. The proportionality 
constant is now k instead of 1/tq of the Bjorken model. We define 

Tf - <^o) (10) 

as a parameter which can be regarded as an equilibration time scale in the same way as tq 
in the scaling case. The initial time t^ s is a parameter which is expected to depend on the 
nuclear size and the collision energy and is adjusted separately to each collision. 

As described above, the model is formulated for symmetric collisions in the center-of- 
mass frame of the collision. Since we are interested in the mass number dependence of the 
initial state and large part of the available data is for asymmetric collisions, we would like 
to modify the parametrization (^) to apply to arbitrary collisions, too. 

In symmetric zero impact parameter collisions the rapidity shift is the same on the target 
and the projectile side and yo(p) = for all values of p. For asymmetric collisions the target 
and the projectile rapidities differ in the overall cm frame of the participating nucleons. 
Also, since we determine the initial conditions from the nuclear thicknesses locally in the 
transverse plane, we have to distinguish between the global cm frame and the cm frame for 
the collision of a row of target nucleons with a row of projectile nucleons at given p. Eq. (0) 
is evaluated in this local cm frame with x y defined using y^ m on the projectile side and yj m 
on the target side 



C p exp(aXy P + bXy P + cx Vi p) + C T exp(— aXy T + bx 2 y T — cx y ^] 
x (1 + x y>T )(l - Xy )P ) , 

where :%t(p) = y/\yl^\ on the target (projectile) side. This ensures the important limit 



that for Ta,b —> oo the system is stopped in the cm frame when we use ([TT]) together 
with (|B|). Note that the distinction of local and global cm frames is necessary also if the 
parametrization is used for nonzero impact parameter collisions of equal nuclei. 

We also have to modify the stopping parametrization to account for the fact that for a 
fixed target size the rapidities of the nucleons of a smaller projectile are, on the average, 
shifted more than those of a bigger projectile and vice versa. This means that we must 
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include both the target and the projectile thickness functions into the parametrization (^j). 
Without this modification the average momentum of projectile nucleons after the collision 
with a target of fixed size would be independent of the projectile size and the total momentum 
of target and projectile nucleons after the collisions would not add up to zero in the overall 
momentum frame. In principle this needs not to be the case since the mesonic degrees of 
freedom can balance the momentum but we do not find it plausible that the center of mass 
of nucleons would appreciable differ from the center of mass of produced mesons. 

It turns out that the momentum balance between the colliding nucleons can be approxi- 
mately satisfied if we replace T4 and T B in Eq. (H) with effective values T| ff and T B S defined 
as 




The factor (Tb/Ta) can be considered as an asymmetry correction for the stopping. The 
choice of a = 0.8 leads to very satisfactory momentum balance within few percent for all 
three investigated asymmetric systems. 

We show in Fig. ^| the initial local baryon density n-Q and the local energy density e(z, p) 
for the five collisions we consider. In Fig. |^a Ub(z, p) is plotted as function of z for p = 0. 
In all cases the stopping leads to approximately Gaussian distribution with the maximum 
at the global center-of-mass position for symmetric cases and slightly shifted towards the 
heavier target side for the asymmetric cases. 

Fig. §d shows the local energy density e(z) for p — 0. It has a Gaussian form resulting 
from Eq.([?D and the linear fluid rapidity profile (||). For asymmetric collisions the maximum 
is at the local yo(p = 0) > = 0. In transverse direction the energy density is proportional 
to Ta(p) Tb(p) resulting in the shapes shown in Fig. |3]c. 

III. RESULTS 

The above-described initial state is evolved in time by solving the ideal fluid hydrody- 
namical equations numerically on a 2+1 dimensional grid described in detail in [|]] where 
different equations of state are compared. Here we show results only for an equation of state 
with phase transition to QGP at T c = 165 MeV labeled as EOS A in [||. 

Freeze-out occurs when the mean free path of particles is of the same order than the size of 
the fireball. For the collisions of light nuclei this criterium leads to freeze-out temperatures 
of Tf pa 140 MeV f9j. Instead of taking contours of constant temperature we define the 
freeze-out on a space-time surface of constant energy density of e f = 0.15 GeV/mi . This 
results in an average freeze-out temperature of Tf w 140 MeV. The freeze-out contours are 
plotted in Fig. [| for the different collisions we have studied. An increase in the size of the 
system leads to longer lifetimes as expected. One sees further that in asymmetric collisions 
the heavier target side lives longer. 
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For Pb+Pb collisions we find indications (see below) that the canonical freeze-out tem- 
perature, Tf ^140 MeV, might be too high for this larger system. Therefore we calculated 
the lead-on-lead collision also with freeze-out at Ef — 0.069 GeV/fm 3 corresponding to an 
average freeze-out temperature of Tf « 120 MeV. 



A. Spectra 

Particle spectra are calculated using the prescription of Cooper and Frye fH| including 



the same hadrons as in the construction of the equation of state f|. In the discussion of 
particle spectra we concentrate on the most abundant particles, i.e. negative particles, pions 
and (net) protons because for them the hydrodynamical treatment with assumed chemical 
equilibrium is most reliable. We show in Fig. |5| the rapidity distributions of negative particles 
for the various collisions. A good agreement is seen in these rapidity spectra. To a small 
extent the agreement in the normalization is improved by adjusting the £ parameter in ([[]) 
which accounts for the impact averaging. This parameter is determined for each collision 
system separately by optimizing the fit to the negative particles, but as shown in Table | 
the variations of £ are quite small. The calculation for Pb+Pb collisions with the lower 
freeze-out temperature of 120 MeV is shown as dotted line in Figs. |]-|8|. 

The rapidity spectra of negative particles are most sensitive to the initial energy dis- 
tribution. Our parametrization (|7]) which includes only two free parameters c £ and a £ in 
the expression of the width (H) reproduces all five measured rapidity spectra of negative 
particles. This suggests that (|7|) could be used for any collision system at SPS energies. 
One should keep in mind, however, that the optimal values for the parameters c £ and a £ in 
Eq. (g) depend on the equation of state used for the hydrodynamical evolution. 

The rapidity spectra of net protons are given in Fig. || For S+S the p — p distribution 
shows two maxima in the target and projectile fragmentation region which are nicely repro- 
duced by the calculation. For Pb+Pb collisions stronger stopping is expected. The choice 
of parameters in Eq. (f|) together with the EOS is able to account for the larger stopping in 
a quantitatively satisfactory way as seen from the comparison with preliminary data from 
NA49 @. 

In the asymmetric collisions the net proton rapidity distributions are reproduced only 
around midrapidity and in the forward region while at target fragmentation region the net 
proton yields are underestimated considerably. The reason for this deviation can be un- 
derstood by considering the geometry of an asymmetric collision. The initial state of the 
hydrodynamical calculation is determined by using the cylindrical volume as cut out of the 
target by the smaller projectile nucleus. However, during the first collisions the projectile 
nucleons moving through the target acquire transverse momentum, some of them penetrate 
into the region outside the cylinder and interact further. Similarly some of the produced 
particles can interact with the nucleons outside the cylinder. Nucleons from these interac- 
tions are not included in the hydrodynamic description but are observed experimentally in 
the target fragmentation region. It is reasonable to assume that these additional protons 
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will not play any role in the hydrodynamical calculation. This assumption is supported by 
the A rapidity distributions in Fig. |6] which show good agreement between the calculated 
and measured spectra even in the target fragmentation region. 

In Fig. [7| we show that the measured transverse momentum distributions of negatives 
and 7r°'s agree quite well with the calculated spectra except in the case of O+Au collisions, 
where the large low-pr enhancement cannot be reproduced and the calculation exhibits 
slightly too much flow. For Pb+Pb the very forward pions are flatter than in the calcula- 
tion. The particle density at this edge of the phase space, however, is becoming so small 
that hydrodynamics with relatively strong longitudinal flow can lead to an artificially large 
transverse cooling. Thus we should not expect that hydrodynamics can describe the far 
edges of fragmentation regions well. 

A closer look on the experimental transverse spectra shows that nearly all of them have 
a small but visible concave shape in my — m which is not reproduced by the calculation. 
The calculation gives only the average slope right. We think that the explanation is the too 
simple freeze-out procedure applied in the calculation. In fact, in such a small system as 
created in the nuclear collisions, the last interaction of a final hadron can take place with 
non-zero probability at any space-time point. Therefore the spectra receive contributions 
from all temperatures and not only from the assumed freeze-out temperature. Note that the 
number of particles in the tail of transverse distribution is very small indicating that rare 
earlier escapes of high-pr particles would not invalidate the hydrodynamic treatment of the 
bulk of the final matter. Improvements have been suggested |19[ but they are difficult to 
implement technically. 

Finally, we show in Fig. ^ the transverse momentum distribution of heavy particles. In 
S+S collisions the net proton distribution is in agreement with the data supporting the 
assumption of freeze-out at a common transverse flow velocity with pions. For S+Ag and 
S+Au we obtain agreement in the py-spectra of A's. The agreement in the normalization 
should be considered accidental. It is obtained with a freeze-out temperature of Tf « 140 
MeV and assuming full chemical equilibration. More detailed studies on chemical equilibra- 
tion indicate higher chemical freeze-out temperatures but not full chemical equilibration of 



strangeness |2D|,2T|. In our calculation the normalization is achieved because the lower freeze- 
out temperature is compensated by the assumption of full chemical equilibration of A's. In 
the transverse spectrum the lower freeze-out temperature is balanced by the transverse flow. 

For Pb+Pb collisions the calculated transverse spectrum of net protons is slightly steeper 
than the data even at central rapidities if a freeze-out temperature Tf ~ 140 MeV is used. 
This was already observed in |22|] and is an indication of the freeze-out taking place at lower 
temperature in the larger Pb+Pb system. Therefore we redid the calculation with a freeze- 
out temperature Tf « 120 MeV as suggested in [p3|-p5|. The result is shown as dotted line 
in Fig. H The agreement for the transverse spectrum of net protons is clearly improved. 
The corresponding m^-spectrum of negatives in Fig. [7] is only slightly modified compared 
to the freeze-out at T « 140 MeV. 

We did a calculation for S+S collisions with Tf = 120 MeV, too. While the pt slope of 
negatives still agrees with the data, that of net-protons, and in particular that of kaons and 
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Lambdas, comes out natter than the data. 

We conclude that our calculations support the kinetic estimates for the freeze-out |§, 
especially its dependence on the size of the system. However, the freeze-out temperatures 
favoured by the p^-spectra are model dependent since we found that results for Pb+Pb 
collisions with T c = 200 MeV and Tf ~ 140 MeV are close to those with T c = 165 MeV and 
T f « 120 MeV. 



B. Parameter Systematics 

We next discuss the A-dependence of various quantities in our model. In Table [I] we show 
parameters which define the collision (upper part) and describe the initial and freeze-out 
state. The parameters £ and Tq S have been adjusted to each collision system. £ is always 
close to one showing that the experiments were triggered to most central collisions and 
include nearly the whole energy and baryon number. The symmetric collisions show a trend 
to a lower value of £ indicating that the impact averaging is more important for these cases. 

The Tq S determines the slope of the initial rapidity profile ( |i~0"D as function of z and thus 
the spatial distributions of densities. A simultaneous fit to the longitudinal and transverse 
spectra determines the optimal Tg ff . It turns out that Tq S is very similar, Tq S ~ 1.3 fm/c, for 
all collisions from S+S to Pb+Pb. However, the sensitivity of final spectra on Tq S is not very 
strong, since the density parametrizations are done in rapidity space. Within a variation of 
±0.2 fm/c in Tq the fits remain acceptable. 

As mentioned already earlier, Tq can be considered to give an estimate of the thermal- 
ization time. On the other hand the initial longitudinal extension can be related to the 
equilibration time for the whole system in the global center-of-mass frame. The longitudinal 
extension of all systems is roughly the same and leads to an estimate of a common initial 
time of around 3 fm/c as counted from the time when the front edges of the nuclei meet 
to the time when all quanta from primary collisions have formed and reached approximate 
equilibrium. The difference compared to t^ s arises mainly from the finite thickness of collid- 
ing nuclei. The numbers are consistent with the findings in the three-fluid hydrodynamical 
studies. They indicate a pressure equilibration time to be of the order of 1.5 fm/c for 



Pb+Pb collisions [|26[] at SPS energies. In the hydrodynamic treatment, carried out in the 
center-of-mass system of participating nuclei, the initial conditions are specified at global 
equilibration time of the whole system which, however, does not enter the calculation as an 
explicit parameter. 

In Table | we give values of the energy density e of the initial fireball formed in different 
collisions. We average over the transverse area in a unit rapidity interval at y = 0, 

Jdppejp) 

£ r — — n ' y L6 ) 

in order to compare with various estimates, in particular with those from the experimental 
transverse energy flow. 
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We expect the A-dependence of the initial energy density to lie between the following two 
crude estimates based on the assumption that e is proportional to the number of sources in 
a nucleus-nucleus collision. An upper limit is provided by the possibility that each nucleon 
in the projectile can interact with any nucleon on its path through the target and v. v. with 



probabilities given, e.g. by the Glauber model [27]. The assumption that the energy of the 
nucleons is degraded independently in each collision with a certain fraction transferred to 
the fireball, leads to the scaling of the energy density as 



e G1 oc T A T B oc R A R B oc (AB) 



1/3 



(14) 



The other extreme is the wounded nucleon model which assumes that each nucleon 
contributes as a single source, see also ||. Then the number of sources on the symmetry 
axis is proportional to T A + T B with a scaling behaviour, which for the mass numbers of 
interest, can be approximated with the form 



_WN 



(xT A + T B oc Ra + Rb 

ocA 1/3 + jB 1/3 = (AjB) 1/6 



(A/B) 1 / 6 + (B/A) 



1/6 



2(AB) 



1/6 



(15) 



It should be kept in mind that the physics of particle production is expected to depend on 
the collision energy. Soft processes can dominate at the SPS energy but at RHIC and in 
particular at LHC the hard and semi-hard production of jets and minijets will probably be 
the main mechanism [^,0 • This means that the mass number scaling of e will be different 
at different energies. 

We investigate the scaling behaviour of the initial energy density by fitting the A- 
dependence with 



e = e (A eS B cS y 



(16) 



The fit to the transversely averaged calculated values Er (open circles) against the product 
of effective mass numbers A eG B eS (cf. Eq. ([[])) is shown in a double logarithmic plot in 
Fig. |. The minimum x -fit to the five collision systems results in 7 = 0.22. 

The result is between the two examples mentioned above (cf. Eq. (II ) and flIB"D), some- 
what closer to the wounded nucleon alternative than the Glauber type behaviour. The 



classical Glauber treatment (14) neglects possible interference effects like the Landau- 



Pomeranchuck-Midgal effect |31| which could lead to a weaker dependence on the mass 
numbers. On the other hand, in the wounded nucleon picture the effects of multiple inter- 
actions may require an increase of the source strength per nucleon leading to faster growth 
than the approximate (AB) 1 / 6 behaviour. 

Our initial energy density may also be compared with the estimates from the experimen- 
tal transverse energy flow (or multiplicity densities) using the Bjorken model for the initial 
volume [HI: 



(dE T /dry) m 



(17) 
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These results for tq = 1 fm/c, taken from [32] and | 3~3[| , are also shown in Fig. ^together 
with the fit, Eq. (|T6|). The x 2 -^ s f° r the er values in our model and the Bjorken estimates 



from experiment are: 



sr = 0.63 GeV/fm 3 (A cff 5 ofr ) ' (18) 
e B j = 0.55 GeV/fm 3 (A cff 5 cfr )°' 16 . (19) 



It is interesting to notice that £Bj has the scaling of (AB) 1 / 6 as expected from the wounded 
nucleon picture. 

We see in Fig. ^ that our average initial energy density is roughly a factor of two higher 
than the Bjorken estimate. This difference reflects the assumed hydrodynamic behaviour. In 
the hydrodynamic expansion the thermal energy is transferred to collective flow motion by 
the expansion work of pressure. When the flow is asymmetric in longitudinal and transverse 
directions the stronger longitudinal flow leads to larger energy transfer into the longitudinal 
degrees of freedom. As a result the energy per unit rapidity is decreased in the central 
rapidity region and increased in the fragmentation regions. If we calculate the initial energy 
density using the Bjorken estimate, Eq. (]T7|) from our calculated final spectra, we obtain 
numbers which are consistent with the estimates from experiment. 

We would like to point out that the scaling of energy density is not very sensitive to the 
parameters used in Eq. (|2J) for the initial baryon distributions. The initial energy density is 
dominantly given by the Tq S and the parametrization of the width a £ (§) in the Gaussian 
energy distribution (|7|). This parametrization is fitted to the rapidity spectra of negative 
particles. Therefore the scaling in the initial energy density is mainly extracted from the 
experimental pion spectra. 

We next discuss quantities which characterize the freeze-out. The lifetime tf of the 
fireball, defined as the freeze-out time at the center, increases as expected with the size of 
the fireball. We see a doubling of the lifetime when going from S+S to Pb+Pb indicating 
approximately linear dependence on the nuclear radius. 

An important quantity from the point of view of the hadron spectra is the transverse flow 
velocity at freeze-out. In a local hydrodynamical calculation the transverse flow is different 
for different fluid cells. In order to compare with average flow velocities of simpler models 
34| or phenomenologically extracted values []35] we have to average over the relevant range 



of the freeze-out surface. At central rapidities we choose for convenience to average over all 
fluid cells with rapidity \y\ < 0.25 in the global fireball rest frame. This corresponds to the 
region of highest transverse flow. The resulting values for (v p ) are given in Table |I[ Note 
that the range \y\ < 0.25 does not necessarily correspond to the rapidity region where the 
transverse momentum spectra in Figs. [7| and |8| are measured. 

The numbers for (v p ) are very similar with the exception of S+S. The fact that we get 
for O+Au a rather similar average flow velocity than for Pb+Pb is a result of taking the 
same freeze-out energy density independent of collision size. Using for Pb+Pb collisions 
the lower freeze-out energy density corresponding to Tf w 120 MeV we find a considerable 
stronger flow than in the smaller systems. A system-size dependent freeze-out can naturally 
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be incorporated if one uses a dynamical freeze-out criteria as suggested in [34| , 24] resulting 
in a decrease of freeze-out temperature with increasing nuclear size. 

The flow values for the symmetric collisions may be compared with the data analysis 
They extracted the mean transverse velocity from a fit to various hadron 
0.24 ± 0.10 for S+S and 



of NA44 35 



in 



spectra obtaining (v p ) = 0.24 ± 0.10 for S+S and (v p ) = 0.36 ± 0.14 for Pb+Pb |35 
agreement with our studies. Our freeze-out temperature Tf ~ 140 MeV agrees with their 
fitted value Tf = 142 ± 5 MeV for S+S. But for Pb+Pb they favour a somewhat higher 
value Tf = 167 ± 13 MeV while our results favour the lower freeze-out temperature of 120 
MeV leading to (v p ) 



0.44 for Pb+Pb in agreement with |23 



IV. CONCLUSIONS 

We have constructed a model for initial conditions in nucleus-nucleus collisions at CERN- 
SPS energies based on the parametrization of baryon stopping by specifying the initial baryon 
number distribution dN^/dy. The parametrization is expressed as a function of the local 
thickness of colliding nuclei in the transverse plane. The major advantages of this approach 
are that it reduces the arbitrariness of choosing the initial conditions for different collisions 
and allows a proper treatment of transverse geometry of the colliding nuclei. We think that 
this is an important improvement compared to the earlier approaches fUH where different 
collisions are not related and the longitudinal dependence is treated equally both on the 
symmetry axis at p = and at the outer edge (p = Rp) of the collision system. 

The success in reproducing the single particle spectra for various collisions gives confi- 
dence that the deduced mass number dependencies of the initial and final state quantities 
are reasonable. An important example is the scaling of the the initial energy density. Even 
though the observed dependence of hadron multiplicity on the mass number is close to linear, 
corresponding to the (AB)^ 1 / 6 ) scaling for the density, the initial energy density grows faster 
as shown in Fig. |9] reflecting the effect of hydrodynamic expansion on the final transverse 
energy flow. For the S+Au and Pb+Pb collisions we obtain a ratio Es+Au/tpb+Pb ~ 2/3 
which, in addition to the volume effects, could be related to the observed features of J/\l/ 
suppression ||36|| . In our parametrization the center of the fireball starts out in the quark- 
gluon plasma phase for all considered collisions. However, we consider the scaling more 
reliable than the absolute value of energy density which is more sensitive on the details of 
the parametrization. One should also note that in smaller systems the high density occurs 
only for a short time interval in a small fraction of the total volume. 

For further improvements in determining the initial conditions of nuclear collisions a 
better understanding of the freeze-out physics is needed, e.g. to allow emission from the 
whole space-time volume |19[ . For a proper treatment also a sequential freeze-out for various 



particle species should be implemented . Since the particle ratios are a result of the whole 
evolution and not only of the freeze-out conditions this would help more directly in pinning 
down the initial state in nuclear collisions. 
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TABLE |. Summary of parameters characterizing different nuclear collisions. A QGP 
equation of state with T c = 165 MeV (EOS A Q) is used. Quantities in the first part 
characterize the colliding system, those in the second part the initial conditions of produced 
matter where values of £ and Tq are obtained from a fit to each specific collision, and the 
quantities in the third part describe the system at the freeze-out. For Pb+Pb collisions the 
numbers in parenthesis correspond to a freeze-out energy density e f = 0.069 GeV/fm 3 with 
an average freeze-out temperature of 120 MeV. The Bjorken estimates EBj are deduced from 
experimental dEx/dr] distributions taken from |32[ and ||33|| . 



collision 


S + S 


O + Au 


S + Ag 


S + Au 


Pb + Pb 


B + A 


32 + 32 


16 + 197 


32 + 108 


32 + 197 


207 + 207 


lab. energy (GeV/A) 


200 


200 


200 


200 


158 


l/cms 


3.03 


2.54 


2.73 


2.62 


2.92 




0.9 


1.0 


0.9 


0.95 


0.9 


r cff (fm/c) 


1.2 


1.4 


1.4 


1.5 


1.3 


e R (GeV/fm 3 ) 


2.5 


4.1 


3.8 


4.8 


6.7 


exp. e B j (GeV/fm 3 ) 


1.3 


2.3 


2.1 


2.6 


3.0 


T R (z = 0) (MeV) 


187 


203 


202 


212 


227 


R(£A) (fm) 


3.15 


6.37 


4.96 


6.37 


6.25 


R{£B) (fm) 


3.15 


2.48 


3.15 


3.28 


6.25 


cent, tf (fm/c) 


6.3 


6.5 


7.5 


8.0 


11.8 (14.7) 


(v p )(\y z \<0.25) (c) 


0.27 


0.34 


0.31 


0.35 


0.33 (0.44) 
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FIGURE CAPTIONS 



Fig. [I|. Rapidity spectra of p — p in Pb+Pb collisions for the Landau initial conditions 
with no initial flow velocity. The calculated spectrum is compared with preliminary data 
from NA49 @. 

Fig. ||. Initial baryon density distribution 7ns in the z-p plane for a S+S collision plotted 
in the overall center-of-mass frame at fixed initial time. 

Fig. [| Initial density distributions for S+S (solid), O+Au (dotted), S+Ag (dot-dashed), 
S+Au (dashed) and Pb+Pb (solid) collisions. In panel (a) the baryon density n B and in (b) 
the energy density e are shown as functions of z on the symmetry axis p = 0. In panel (c) 
the radial dependence of energy density e is shown at z — 0. 

Fig. Freeze-out contours £/ = 0.15 GeV/fm 3 for the various collisions. The lines are as 
in Fig. H 

Fig. H Comparison of rapidity spectra of negative particles measured in various collisions 
with our results for EOS A @. The average freeze-out temperature is 140 MeV (solid line) 
and 120 MeV (dotted line). The data sets are: S+S (NA35) Q; Pb+Pb (NA49) |T2j; O+Au 
(NA35) 0; S+Ag (NA35) and S+Au (NA35) 0. The Pb+Pb data is preliminary. 

Fig. p. Comparison of rapidity spectra of p — p and A measured in various collisions with 
our results for EOS A |4]]. The average freeze-out temperature is 140 MeV (solid line) and 
120 MeV (dotted line). The data sets are: p-p in S+S (NA35) fTPf; p-pin Pb+Pb (NA49) 
@; p - p in S+Ag (NA35) p - p in S+Au (NA35: y > 3) (T§ (two data points with 
y < 3 are proton data from NA44 pfl); A in S+Ag (NA35) |L7| and A in S+Au (NA35) 



T7| . The Pb+Pb data is preliminary. 



Fig. |7]. Comparison of transverse momentum spectra of negative particles and ir° measured 
in various collisions with our results for EOS A 0. The average freeze-out temperature is 140 
MeV (solid line) and 120 MeV (dotted line). The data sets are: S+S (NA35) Q; Pb+Pb 
(NA49) |]12| for rapidity intervals of width 0.5 and center at (top to bottom) 3.4, 3.9, 4.4, 
4.9, 5.4 and data sets successively scaled down by 10~ n ; O+Au (NA35) @; S+Ag (NA35) 
for rapidity intervals (top to bottom) 0.8 < y < 2.0, 2.0 < y < 3.0, 3.0 < y < 4.0, and 
4.0 < y < 4.4 and data sets successively scaled by 10~ n (n = 0,1,2,3); S+Au (WA80) @. 
The Pb+Pb data is preliminary. 

Fig. Comparison of transverse momentum spectra of p — p and A measured in various 
collisions with our results for EOS A B. The average freeze-out temperature is 140 MeV 
(solid line) and 120 MeV (dotted line). The data sets are: p — p in S+S (NA35) [0; p — p 
in Pb+Pb (NA49) |12j for rapidity intervals of width 0.5 and center at (top to bottom) 2.9, 
3.4, 3.9, 4.4, 4.9, 5.4 and data sets successively scaled down by 10~ n ; A in S+Ag (NA35) 
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17j and A in S+Au (NA35) |L7|]. The Pb+Pb data is preliminary. 



Fig. ^. Initial average energy density plotted against A cR B cS . Open circles show the cal- 
culated values Eq. (|T3[) and the diamonds the Bjorken estimates Eq. (|T7|) from experimental 
data p2| , p3| with lines corresponding to minimum \ 2 fits of the form Eq. (|16|). 
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